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We present a comparison between space charge calculations and direct measurements of the trans- 
verse phase space for space charge dominated electron bunches after a high voltage photoemission 
DC gun followed by an emittance compensation solenoid magnet. The measurements were performed 
using a double-slit setup for a set of parameters such as charge per bunch and the solenoid current. 
The data is compared with detailed simulations using 3D space charge codes GPT and Parmela3D 
with initial particle distributions created from the measured transverse and temporal laser profiles. 
Beam brightness as a function of beam fraction is calculated for the measured phase space maps 
and found to approach the theoretical maximum set by the thermal energy and accelerating field at 
the photocathode. 

PACS numbers: 29.25.Bx, 29.20.Ej, 29.27.Bd, 52.59.Sa, 52.65. Cc, 52.65. Rr, 52.70.Nc 



I. INTRODUCTION 

The generation of high-brightness electron beams re- 
mains the principle challenge for a number of linear ac- 
celerator based projects including the Energy Recovery 
Linac (ERL). Design of the electron sources relies heavily 
on the use of space charge simulations P, [H, H, 0, H, [if- 
Still, there remain a number of open questions with re- 
gards to understanding and modeling of space charge 
dominated bunched beams in photoemission guns. Dif- 
ferent space charge codes use varying degree of approxi- 
mations to capture the most significant physics relevant 
to beam dynamics in photoinjectors. Certain codes allow 
self-consistent inclusion of complex conducting bound- 
aries at the expense of considerable increase in the re- 
quired computation time 0, [H , while many of the main- 
stream codes widely used in the design of photoinjectors 
assume open boundary conditions everywhere except in 
the vicinity of the photocathode. Due to an inability to 
use direct self-force calculations in a bunch with ~ 10 9 
particles, artificial smoothing of the space charge poten- 
tial is employed when tracking macroparticles, which rep- 
resent the actual bunch, either through meshing of the 
electron cloud or via introduction of an effective size to 
the macroparticles in a point-to-point calculation. As 
a result, and depending on the simulation parameters, 
the space charge force in simulations can either be overly 
smooth or grainy as compared to the actual self-force in 
space charge dominated beams. Additional assumptions 
are employed by different codes with respect to details 
of modeling the emission process from the photocathode, 
whether or not the velocities of individual electrons in 
the rest frame of the bunch are treated as negligible, 3D 
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or 2D or lower dimensionality nature (e.g. uniform cylin- 
ders of HOMDYN) of the space charge, etc. The validity 
of these assumptions must be evaluated for each individ- 
ual case. Yet, data comparing directly measured beam 
conditions, detailed phase space distributions in particu- 
lar, with simulations is sparse for space charge dominated 
bunched beams such as found in either DC or RF pho- 
toinjectors [1, EH [HI] . Even a relatively simple configu- 
ration involves a number of "free" parameters that need 
to be varied within the uncertainty of the measurement 
such as the RF phase in an RF gun in order to obtain 
good agreement between simulations and measurements 

We present direct measurements of the transverse 
phase space distribution using a simple beamline con- 
sisting of a DC gun followed by an emittance compen- 
sation solenoid. The number of potential variables af- 
fecting beam performance is reduced to the bare min- 
imum in such a setup. Careful characterization of the 
initial conditions such as the laser transverse and tempo- 
ral profiles and thermal emittance of the photocathode 
allows us to carry out a direct cross-checking between the 
measurements and 3D space charge simulations using the 
codes Parmela3D [l[ and GPT [J. Phase space distribu- 
tions are then used to calculate beam brightness vs. the 
included beam fraction, which is compared to the the- 
oretical limit set by the thermal transverse energy and 
accelerating field at the photocathode. 

The paper is organized as following: Section |TT] de- 
tails the experimental setup and beam diagnostics used 
in this work as well as our experimental procedures. Sec- 
tion IIIII presents simulation details and data processing 
procedures used to extract information from the mea- 
sured data such as second moments and rms emittances. 
Comparison between data and simulations follows in Sec- 
tion [TV] Finally, we conclude with a discussion and out- 
look for future work. 
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FIG. 1: Beamline used in the space charge studies. Beam direction is to the left. 



II. EXPERIMENTAL SETUP 

A. Beamline 

Fig. Q] shows the experimental setup, which consists of 
a high-voltage DC gun [12j followed by a solenoid used 
for emittance compensation located at 0.335 m from the 
photocathode to its center. Fig. [2b shows the distribu- 
tion of the magnetic field in the solenoid as calculated 
by POISSON [13[ and the actual measured values us- 
ing a Hall-probe. Both the emittance measurement sys- 
tem (EMS) and an insertable viewscreen are positioned 
1.244 m from the photocathode. Additionally, the beam- 
line is equipped with a deflecting RF cavity [3| used to 
characterize the initial temporal profile of photoemitted 
electrons, two beam scanners and a Faraday cup all in- 
tegrated into a data acquisition system for direct phase 
space measurement. Details on the EMS are in the next 
subsection. Two different types of materials have been 
employed for viewscreens: high sensitivity BeO used in 
temporal measurements with very low bunch charges and 
the deflecting cavity and CVD diamond used with aver- 
age beam currents of up to 100 fj,A. Each viewport is 
equipped with 12-bit CCD camera interfaced to the con- 
trol system. 
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FIG. 2: Field profiles for (a) the DC gun at 250 kV and (b) 
the solenoid with 5 A excitation current. 



The HV DC gun, initially designed for 750 kV, was 
operated at 250 kV. The field distribution is shown in 
Fig. [2^. The gun has reached 420 kV voltage during 
high voltage processing in the year 2007. Since then, 
however, we had to limit the gun voltage to a conservative 
value below 300 kV due to field emission problems. Upon 
disassembly of the gun we found a considerable amount 
of dust coming from the ceramic resistive coating, which 
is believed to be the primary reason behind the strong 
field emission. Work is underway to eliminate this source 
of dust from the gun. 

The laser system has been detailed elsewhere [l5|, [l6[ . 
The laser spot-size was monitored using a 12-bit CCD 
camera at the location of a virtual cathode. A Pockels cell 
was used to reduce the 50 MHz train of pulses, with an 
average power of about one Watt and 520 nm wavelength, 
to a lower duty factor for beam measurements. A typical 
average current during phase space measurements was 
between 10 to 100 /iA. 



B. Emittance Measurement System 

Considerable care is required when designing a sys- 
tem for direct measurement of the phase space of space 
charge dominated beam [U HH • The EMS employed in 
this study is a double-slit system, with non-moving parts; 
see Fig. [3J The beam motion is achieved with a pair of 
corrector coils, designed with a vanishing sextupole com- 
ponent in order to provide uniform kick across most of 
the vacuum beampipe cross-section (see Fig. [4j. Each of 
the two scanners, one prior to the 1st slit, and the other 
between the two slits consists of a pair of identical coils 
with opposite direction of excitation current. The coils 
in each beam scanner have been measured to cancel each 
other to better than 1%. Thus, each scanner changes 
only the position of the beam without affecting its di- 
vergence. Additionally, a pair of horizontal and vertical 
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FIG. 3: (a) Emittance Measurement System. One beam scan- 
ner (not shown) precedes the 1st slit, (b) Details of the 1st slit 
showing water-cooled 200 fim armor slit and 20 /jm precision 
slit. 



day cup configurations. Monte-Carlo analysis on scat- 
tered radiation using GEANT4 [19( has been carried out 
for the full EMS system showing excellent signal to noise 
performance of the system with cither of the two config- 
urations [20l | . 

All of the measurements reported in this work have 
been done using a double-slit method. The maximum 
scan rate for the beam scanners was 200 Hz. The signal 
from both slits has been detected using the Faraday cup 
connected to a low noise current amplifier. A solenoid 
positioned just before the Faraday cup was used to focus 
the beamlet to the center of the Faraday cup for increased 
charge collection efficiency. A typical transverse phase 
space scan of 100 x 100 points would take on the order of 
one minute. 

The slit opening size and distance between the two slits 
have been determined by solving coupled beam envelope 
equations for a beamlet selected by the 1st slit with in- 
clusion of the space charge force: 
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Here / is the beam peak current after passing through 
the 1st slit, Iq — 17 kA is the Alfven current, and ((3j) 
is the normalized momentum. After a vertically select- 
ing slit with opening d small compared to the beam 
size, one has a y = d/\/Y2 and normalized rms emit- 
tance e n<y = e rit yo(d/y/ll2))/a y o, where e n ,„o and a yQ 
are the emittance and vertical size of the full beam prior 
to the slit. By solving Eq. Q]for the beamlet size u y (L) 



FIG. 4: Uniformity of B-field integral for beam scanner 
coils. Quantity being plotted is T(x, 0)/X(0, 0) — 1 and 
T(0, y)/l(0, 0) — 1, where T(x,y) = J B x (x,y, z)dz and x — 
y = corresponds to the center of the beampipe. 



steering coils before each of the slits allows correction for 
yaw/pitch alignment errors, leading to overall relaxed tol- 
erances for the device. Finally, a small solenoid (about 
±1° rotation angle) is positioned between the two slits 
to allow for roll compensation, although its use proved 
unnecessary in practice. 

The 1st EMS slit consists of an armor slit with a 200 
fim opening followed by a 20 /im vertically selecting preci- 
sion slit brazed to the water cooled armor slit; see Fig.[3p. 
With this design most of the beam power is intercepted 
by the armor slit. ANSYS analysis shows that the EMS 
system is able to perform without significant deformation 
of the precision slit (< 10%) with 1 kW of incident beam 
power. 

Phase space measurements can be carried out using 
either single slit and viewscreen or double-slit and Fara- 
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FIG. 5: Emittance overestimation due to space charge for 
different slit openings. Beam parameters: tn,x, y = 0.3 /xm, 
rms bunch duration a t = 20 ps, a X]V — 1.3 mm, total rms 
divergence PT,x,y = 0.72 mrad, kinetic energy 0.5 MeV, and 
bunch charge 80 pC. Red dot shows the actual separation 
between the two slits (38 cm). 
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FIG. 6: Typical transverse laser profile. Black solid lines show 
projected horizontal and vertical profiles. 
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FIG. 7: Temporal profile of electron bunches as measured by 
the deflecting cavity with negligible charge per bunch. Red 
solid line shows the reconstructed profile used in simulations. 
Refer to the text for details. 



at the location L of the 2nd slit (or the viewscreen) , and 
comparing a y (L)/L to the uncorrelated divergence a' yQ 
at the 1st slit position, one can gauge the effectiveness 
of the slit system in the presence of space charge. Here 
we note that the simple requirement that the beamlet 
after the 1st slit be emittance dominated [TtJ is neces- 
sary but not sufficient. Because the ratio of the space 
charge and emittance terms in Eq. (1) scales as fast as 
tx <jy, meeting such condition at the location of the 1st 
slit does not ensure that the beamlet stays emittance 
dominated all the way through to the location of the 2nd 
slit. Numerically solving the coupled equations Eq. [JJ 
on the other hand, allows proper characterization of the 
slit performance. Fig. [5] shows the results of emittance 
overestimation by the double-slit method in the case of 
80 pC bunches with the parameters as indicated in the 
figure. The separation between the two slits was chosen 
to be 38 cm leading to the overestimation in emittance 
measurement to be less than 10% of 0.3 fim normalized 
rms emittance at 0.5 MeV kinetic energy. 

C. Experimental Procedures 

Measurements have been taken at 3 different bunch 
charges: 80 pC, 20 pC and 0.5 pC. The measured laser 
intensity stability was 2% rms. The laser spot was ini- 
tially magnified, then passed through a 2.6 mm diameter 
aperture, which was 1:1 imaged onto the photocathodc. 
The laser pointing stability was 60 /im rms in each trans- 
verse direction. Each data set involved taking multiple 
images of the laser spot on the virtual photocathode, and 
an image representing the average centroid position was 
chosen for simulations as detailed in the next section. A 
typical transverse laser spot profile is shown in Fig. [6l 

The laser pulses were temporally stacked using three 



birefringent crystals [16j |. Direct measurement of the 
initial temporal distribution of electrons was performed 
with the deflecting cavity with negligible charge per 
bunch. See Fig. [7J The resolution of the temporal mea- 
surement in this case is 1.5 ps rms as limited by the RF 
to laser synchronization and finite electron beam spot 
size. To obtain the actual temporal profile of the electron 
distribution, the data was fitted with 8 Gaussians, then 
each Gaussian was assigned 1.0 ps sigma corresponding 
to the value found in an autocorrelation measurement 
for an unshaped laser pulse [l5[ • Both the fit to the data 
and reconstructed temporal profile used in simulations 
are shown in Fig. [JJ 



Early on in the measurements, we were able to observe 
asymmetric transverse phase space distribution, e.g. see 
Fig. [TS1 To eliminate possible causes for such an oc- 
currence, careful beam based alignment was carried out 
before each data collection. A small laser aperture (0.25 
mm) placed concentrically with the larger one was used 
to create a small beam with negligible charge per bunch. 
The beam centroid vs. the solenoid current data was fit 
to obtain both the angle and the offset of the magnet's 
magnetic axis with respect to the beam. The solenoid has 
been physically adjusted so that its axis coincides with 
the beam axis to within a few 10's /im and /trad in offset 
and angle respectively. Similarly, the center of the aper- 
ture was imaged to coincide with the electrostatic center 
of the gun to about 10 /im, ensuring that the central orbit 
is well aligned throughout the system. 
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FIG. 8: Comparison of beam envelope vs. longitudinal posi- 
tion calculations for Parmela3D (solid lines) and GPT (open 
symbols) for 80 pC charge per bunch. Dashed line shows 
Parmela calculations using 2D space charge routine. The 
solenoid current is 3.6A. The gun voltage is 250 kV. 
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FIG. 9: Comparison of emittance vs. longitudinal position 
calculations for Parmela3D (solid lines) and GPT (open sym- 
bols) for 80 pC charge per bunch. Dashed line shows Parmela 
calculations using 2D space charge routine. The solenoid cur- 
rent is 3.6A. The gun voltage is 250 kV. 



III. SIMULATION PARAMETERS AND DATA 
PROCESSING PROCEDURES 

A. Simulation Parameters And Conditions 

As the quantum efficiency (QE) of the photocathode in 
the region of interest was found to have about 10% peak- 
to-peak fluctuations (in absolute terms the QE for GaAs 
was about 6% in this work), the measured laser trans- 
verse profiles and the temporal shape shown in Fig. [7] 
were used to created 3D distributions for simulations. 
The thermal emittance of GaAs photocathodes has been 
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FIG. 10: Comparison of emittance vs. position calculations 
for GPT in case of different boundary conditions: open (06), 
approximate (A6), and Dirichlet (D6). Refer to the text for 
details. Bunch charge is 80 pC, the solenoid current is 3.6A. 
The gun voltage is 250 kV. 100k macroparticles were used in 
the simulations. 



investigated previously [21| and for an rms laser spot 
size <7_l was found to be e r . 



<j± y/ kT± I mc 2 with 
kT± = 120 ± 8 meV for 520 nm light, mc 2 is electron 
rest energy. 

We have used two 3D space charge codes: Parmela3D 
which employs a fast Fourier Transform method for solv- 
ing the Poisson equation on a 3D grid [12] and Gen- 
eral Particle Tracer (GPT) with a non-equidistant mesh 
solver for the space charge force calculation 23]. The 
same field maps for both the gun and the solenoid magnet 
were used in both codes. Convergence of the calculation 
results has been checked for 20k, 100k, and 500k par- 
ticle distributions and different mesh sizes. Results for 
the beam envelope are presented in Fig. [HI Fig. [5] shows 
the results for emittance vs. longitudinal position for the 
case of 80 pC bunches. An additional difference between 
Parmela and GPT is that the former reports the rele- 
vant beam parameters as a function of time, whereas 3D 
coordinates of the bunch were projected to a given longi- 
tudinal position for the latter. We observe that sufficient 
convergence is demonstrated with 100k macroparticles 
for both Parmela3D and GPT. The mesh size was set to 
64 x 64 x 64 for Parmela3D and 50 x 50 x 50 for the non- 
equidistant mesh Poisson equation solver in GPT with 
the bounding box size set to 5a in each dimension. Ad- 
ditionally, the Poisson equation solver in GPT provides 
a choice of 3 different boundary conditions at the bound- 
ing box: Dirichlet with zero potential, an open boundary, 
and an approximate boundary in which the potential at 
the bounding box is assigned analytically computed val- 
ues from a uniform elliptical cylinder with rms dimen- 
sions set equal to those of the actual bunch. (Zero po- 
tential is assigned at the cathode to include image charge 
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FIG. 11: Example of measured transverse phase space (a), 
with a contour map of a binary image after applying 0.6% (b) 
and 4.4% (c) threshold of the peak intensity. 
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FIG. 12: Noise subtraction verification procedure, (a) Ex- 
ample transverse phase space with the contour (1) obtained 
through the boundary detection algorithm and the grown 
contour (2), which corresponds to 50% of the available data 
treated as noise, (b) Normalized rms emittance calculation as 
a function of included area after the noise subtraction. Refer 
to the text for details. 



effects for all boundary condition types). Comparison be- 
tween different boundary conditions is shown in Fig. [TU] 
for emittance vs. longitudinal position. Beam envelope 
calculations for these boundary conditions produce essen- 
tially identical results for these parameters, and therefore 
are not shown. Open boundary conditions were used for 
subsequent comparisons with the measurements. 



B. Data Processing Procedures 

In order to extract second moments from the measured 
beam profiles and phase space distributions, the data re- 
quires appropriate noise subtraction. The general ap- 
proach follows the general notion of self-consistent unbi- 
ased rms emittance analysis (SCUBEEx) (24j: i) a certain 
contour delineates the signal plus noise region from the 
noise only region; ii) the average intensity of the outside 
region represents a noise bias; iii) the noise bias is sub- 
tracted from the data while the outside region is assigned 
intensity; iv) the contour is grown, and the parameter 
of interest should not change significantly as a function 
of the included area once all signal is accounted for if uni- 
form random noise is present. We have used two types 
of contours: circular type for the viewscreen data, and 
a special boundary detection technique for the measured 
phase-space distributions. The boundary detection tech- 
nique is based on the following observation: a binary im- 



SOL1 = 3.7A 

FIG. 13: Comparison of measured and simulated vertical 
transverse phase space distributions for 0.5 pC bunches at 
z = 1.244 m. Solenoid current is 3.7 A. Corresponding rms 
normalized emittances: e n ,y ~ 0.31 ± 0.04 /jm (data), 0.29 
fim (Parmela3D), 0.28 (im (GPT). Corresponding rms sizes: 
cr y = 1.15 ± 0.05 mm (data), 1.14 mm (Parmela3D), 1.14 
(GPT). 



age obtained by applying a threshold to the phase space 
2D distribution is likely to form a continuous region for 
the signal, and many individual islands for the noise. See 
Fig. [TTJ The boundary detection algorithm proceeds as 
follows: i) the data is convolved with a n x n square 
(image blurring); ii) the smallest threshold is found that 
generates a single continuous island; iii) n is incremented 
and step i) is repeated. The process stops when the island 
starts to include chunks of noise region, which becomes 
clearly visible. Once the boundary has been found, noise 
subtraction is verified by growing/shrinking the contour. 
The contour growing is stopped when less than half of the 
whole image area becomes available for noise estimation. 
Fig. [T2] illustrates the procedure further. The change in 
the parameter of interested (e.g. rms emittance) vs. the 
included area delineating the signal from noise regions 
represents the uncertainty in the measurement due to the 
noise subtraction. Emittance or rms values so calculated 
correspond to 100% of the beam. 



IV. COMPARISON OF MEASUREMENTS 
WITH SIMULATIONS 

A. Transverse Phase Space Distributions 

Fig. [13] shows the comparison of measured transverse 
phase space at the location of the 1st slit (z = 1.244 m) 
for 0.5 pC charger per bunch with Parmela3D and GPT 
simulations. The resolution of the measured transverse 
phase space is 90 x 90 steps. In the case of simulations, 
each image is produced using a 300 x 300 2D histogram 
with additional convolution (blurring) with a 3x3 square. 
Each image is normalized to the same maximum inten- 
sity value. An identical color map to that of Fig. [6] is 
used throughout. As expected, the calculated rms nor- 
malized emittance is in good agreement with the thermal 
emittance value for this case. 

Fig. [14] and [15] show the comparison for 20 pC and 
80 pC charge per bunch respectively vs. solenoid lens 
strength. The streak features seen in the measured phase 
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FIG. 14: Comparison of measured and simulated vertical transverse phase space distributions for 20 pC bunches at z = 1.244 m. 
Data representing measurements, Parmela3D and GPT calculations is arranged in rows with different strength of the solenoid 
lens corresponding to column position. 
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FIG. 15: Comparison of measured and simulated vertical transverse phase space distributions for 80 pC bunches at z = 1.244 m. 
Data representing measurements, Parmela3D and GPT calculations is arranged in rows with different strength of the solenoid 
lens corresponding to column position. 
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FIG. 16: Comparison of vertical rms beam size (a) and nor- 
malized rms emittance (b) at the location of the measurement 
z = 1.244 m vs. solenoid excitation current for 20 pC bunches. 
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FIG. 17: Comparison of vertical rms beam size (a) and nor- 
malized rms emittance (b) at the location of the measurement 
z = 1.244 m vs. solenoid excitation current for 80 pC bunches. 



space are due to the motion of the laser spot at the pho- 
tocathode. Good qualitative agreement can be seen for 
20 pC/bunch data, while some discrepancy in the shape 
of transverse phase space distributions can be seen at 
larger solenoid current values for 80 pC/bunch. 




FIG. 18: Example of asymmetric phase space distribution 
measured for 80 pC/bunch at the solenoid current of 3.7 A 
(left) and 3.8 A (right). The top part of the asymmetric tail 
was truncated during the measurement. 
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FIG. 19: Calculated vertical transverse phase space using 
GPT (a) at the location of the 1st slit (z = 1.244 m) for 
a given transverse initial laser spot size (b). Transverse pro- 
file is shown at the location inside the solenoid (c), z = 0.35 
m. Note two different scales for (b) and (c). The portion of 
the distribution responsible for the tail is also shown. 



C. Asymmetric Phase Space Distributions 



B. Second Moments Of The Beam 

Fig. [16] shows a comparison of vertical rms beam size 

(a) at the location of the viewscreen (z = 1.244 m) and 
rms normalized vertical emittance for 100% of the beam 

(b) as a function of solenoid current for 20 pC bunches. 
Excellent agreement is seen for the spot size compari- 
son and good overall agreement for the emittance values, 
although the measured rms emittance appears to be sys- 
tematically smaller for the 20 pC/bunch case. 

Fig. [T"7l shows similar results for 80 pC charge bunches. 
Different sets of curves for EMS and viewscreen simu- 
lations correspond to different laser spots as registered 
for the two data sets. Good agreement between simula- 
tions and measurements of the beam size is seen for 80 
pC/bunch before the formation of a beam waist at the 
location of measurement, while the agreement at larger 
solenoid currents appears to be less conclusive. 



Asymmetric phase space distributions have been mea- 
sured on multiple occasions for space charge dominated 
beam conditions in our setup. For example, see Fig. 1181 
No such asymmetry was observed for low bunch charge 
running under otherwise identical operation conditions. 
It was important to understand the origin of this behavior 
for its subsequent mitigation. We were able to reproduce 
similar phase space distributions in 3D simulations for 
laser spots with noticeable asymmetry. To elucidate the 
mechanism for this tail formation, the particles compris- 
ing the tail have been tagged and their portion of the 
distribution is shown in Fig. [TO] for the transverse profile 
at the location of the photocathode, z = (b), and inside 
the solenoid, z = 0.35 m (c). The asymmetry in the laser 
spot causes the space charge forces to push these par- 
ticles away from the central axis (three is about 2 mm 
difference between the top and the bottom edges of the 
transverse distribution in Fig. [TOb ) so that the particles 
experience a stronger focusing kick from the solenoid lens. 
These particles then undergo a cross-over and form the 
observed phase space tail. Improving the transverse laser 
shape would reduce the asymmetry in the phase space. 
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core e = 0.19 urn, beam fraction = 71% 



core e = 0.18 urn, beam fraction = 66% 
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FIG. 20: Normalized rms emittance vs. included beam frac- 
tion for measured (a) and calculated by GPT (b) phase space 
distributions for 20 pC/bunch. A corresponding Gaussian 
beam is also shown for comparison. The 100% normalized 
rms emittance is 0.43 ±0.05 /im for the measurement (a) and 
0.49 /im for GPT (b). 



core e = 0.44 urn. core fraction = 




FIG. 21: Normalized rms emittance vs. included beam frac- 
tion for measured (a) and calculated by GPT (b) phase space 
distributions for 80 pC/bunch. A corresponding Gaussian 
beam is also shown for comparison. The 100% normalized 
rms emittance is 1.8 ± 0.2 /im for the measurement (a) and 
1.8 /im for GPT (b). 



D. Theoretical Limit To Beam Brightness 

It is instructive to consider rms emittance as a func- 
tion of the contributing beam fraction [l8|. For exam- 
ple, Fig. [20] shows rms normalized emittance vs. beam 
fraction as measured and simulated using GPT for 20 
pC/bunch for a solenoid current of 3.8 A. Fig. |2"T1 shows 
corresponding results for 80 pC /bunch for a solenoid cur- 
rent of 3.6 A. Dot-dashed lines show the expected curve 
for a Gaussian distribution in the phase space with the 
same rms emittance as 100% of the actual beam. It 
can be seen that in the case of 80 pC/bunch, the beam 
strongly deviates from Gaussian distribution, having a 
substantially brighter core. Core emittance, defined as 
£ n , y ,core = de n , y (t; = 0)/c!£ with e niJ/ (£) being the nor- 
malized rms emittance as a function of beam fraction 
< £ < 1, is given in both figures along with the 
core fraction £ core defined as the fraction of the beam 
with the emittance equal to the core emittance value: 

€n,y(,£,core) — ^n,y,core- 

The beam brightness available from photoinjectors 
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FIG. 22: Beam brightness normalized per single bunch as a 
function of contributing beam fraction for 80 pC charge per 
bunch. Corresponding result for an equivalent Gaussian beam 
with the same 100% rms emittance is shown for comparison 
(dot-dashed line). Dotted line shows theoretical limit as given 
by Eq. [5] The rms normalized emittance is 1.8 ± 0.2 fj,m. 



forms through an interplay of several phenomena such 
as space charge dominated beam dynamics in the pres- 
ence of time transient and position dependent external 
fields. The upper limit, however, is set by the thermal 
omittance of the photocathode and the available acceler- 
ating gradient. Consider a short laser pulse illuminating 
a photocathode placed in the accelerating field E cat h ■ 
The electron bunch after the emission will assume a pan- 
cake shape provided that the laser pulse duration is suffi- 
ciently short: a t <C \J <jj_m/eE cat h, rn and e are electron 
mass and charge respectively. This condition is satisfied 
in most operating photoinjectors with bunched beams. 
The maximum charge density that can be supported by 
the electric field is then given by 



dq - , P 

— — — eQ-Cjcath- 

dA 



(2) 



Note the inclusion of the image charge. The average (nor- 
malized) beam brightness can be defined as a ratio of 
average current I avg over its 4D-volume A4 defined for 
{x,p x /mc,y,p y /mc) coordinates: 



B — 



Aa ' 



(3) 



This quantity would be related to x-ray brightness for a 
properly matched undulator, for example. Beam bright- 
ness normalized per single bunch is given by B n ^ avg / f = 
q/Ai, with / being the repetition rate, and q the charge 
contained in the 4D-volume Aa- E.g. the 4D volume 
A4 — dxdp x dydp y / (mc) 2 for a 4D-hypercuboid element 
with sides dx, dp x , dy 7 and dp y . The charge contained in 
A4 can be written as 



—t- dxdp x dydpy — !— . 

kct 2 , 



dA 



(4) 
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with a p being the rms value of the transverse momentum 
(assumed to be isotropic for both transverse directions), 
which is <T. p = \fmkT±_ for Maxwcll-Boltzmann distribu- 
tion of velocities. The dimensionlcss coefficient k depends 
on details of momentum distribution, e.g. k — Ait cor- 
responds to a uniform circular distribution in p x and p y 
with a diameter 4<t p , while k = 2-7T corresponds to the 
peak of a 2D Gaussian distribution. Combining Eq. 0] 
with the charge density as given by Eq. [2l we find the 
maximum beam brightness normalized per single bunch 

Bn.avg _ £o m C 2 E cat h 

This result shows the maximum beam brightness avail- 
able from a photoinjector to be independent from the 
bunch charge, and is determined by the accelerating field 
Ecath and transverse thermal energy kT±_ of the electrons 
leaving the photocathode. To compare Eq. [5] with the 
measured data, we compute the beam brightness per sin- 
gle bunch as 

/ meas. \^n,y(0 J 

where is 100% charge, and e n ,j/(0 is emittance vs. beam 
fraction curve (cf. Fig. [2TJ Eq. refeq:brmeas addition- 
ally assumes an axially symmetric beam with uniform 
phase space distribution inside an equivalent ellipse with 
4Tre n ,y(0 area. Fig. 1^21 shows brightness normalized per 
single bunch vs. the beam fraction. The theoretical (av- 
erage) brightness maximum as given by Eq.[5]is shown as 
well (k = Air). Additionally, {B) n>avg / f is computed for 
a beam that has the same 100% rms emittance as the ac- 
tual beam but adopts a 2D Gaussian distribution in the 
phase space. It is seen that an equivalent Gaussian beam 
does a poor job of describing the 80 pC/bunch beam as 
the core is substantially brighter for the measured beam 
approaching the limit given by Eq. [5l 

V. DISCUSSION AND OUTLOOK 

Benchmarking of 3D space charge codes has been 
performed with the direct measurements of the trans- 
verse phase space for a bunched beam in the space- 
charge dominated regime from a DC gun. Overall, good 
agreement has been found between the measurements 
and simulations. We observe that an equivalent Gaus- 
sian beam assigned the measured 100% rms emittance 
poorly describes the peak brightness available in the 
beam at 80 pC bunch charge due to the presence of a 
substantially brighter core. In particular, for the case 
of the minimum measured 100% rms normalized emit- 
tance € rhy — 1.8 ± 0.2 /im, the core emittance is found 
to be e n ,y,core = 0.31 ± 0.04 /im with 60% beam fraction 
contained in the core. Additionally, the comparison of 
the measured beam brightness vs. beam fraction shows 
that it approaches the maximum theoretical brightness 




0.3 1 , , i 

200 400 600 800 

gun voltage (kV) 

FIG. 23: Calculated normalized rms emittance (100% of the 
beam) at z = 1.244 m vs. the gun voltage for the beamline 
setup used in this study. Refer to the text for details on 
simulation parameters. 



as set by the available accelerating gradient and trans- 
verse thermal energy of the photocathode. While the 
peak brightness of the beam cannot be improved without 
changing kT± and E cat h, it should be possible to bring 
a larger portion of the electron beam to approach the 
brightness limit given by Eq. [5] through proper control of 
the space charge forces. For example, Fig. |2"31 shows sim- 
ulated possible rms normalized emittance (100% of the 
beam) from the same beamline as used in this experiment 
for 80 pC bunches as a function of gun voltage. A uni- 
form cylindrical laser distribution has been used in these 
calculations, with the (same) thermal transverse energy 
corresponding to GaAs illuminated by 520 nm. The laser 
pulse duration in these simulations was 12 ps rms, about 
50% longer that what was used in the measurements re- 
ported here. In addition to employing a longer laser pulse 
and continuing the work to reach the gun design volt- 
age of 750 kV, noticeable improvements are sought for 
the transverse laser shape as well as the pointing sta- 
bility. E.g. based on simulations, the worst 10% emit- 
tance beam fraction in our measurements can be mostly 
attributed to the less than ideal transverse laser profile 
(see Fig. El). 

Finally, we note that despite a significantly more com- 
plicated setup and beam dynamics in the full ERL injec- 
tor [25| , where the bunch undergoes acceleration to over 
10 MeV, bunch compression and matching into the linac, 
a much simpler beamline such as the one considered in 
this study allows exploration of the best beam brightness 
achievable from the whole photoinjector. In particular, 
simulations for the full injector [26|], where substantially 
shorter bunches are produced though subsequent drift 
bunching (~ 3 ps rms) indicate rms normalized emit- 
tances at ~ 1 1 MeV which are about 50% lower but oth- 
erwise very similar to that shown in Fig. 1231 Thus, to 
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continue the work on improving the HV DC gun design 
and pushing for lower omittances, it is sufficient in many 
ways to have a simpler setup with beam diagnostics ded- 
icated to such research similar to the one described in 
this work. 



VI. ACKNOWLEDGMENTS 

We acknowledge control system support by John Dob- 
bins, John Barley, and Mike Forster, laser support by 



Dimitre Ouzounov and Heng Li, magnetic measurements 
of the solenoid by Alexander Temhykh. The following 
individuals are acknowledged for many useful discussions 
and their interest in this work: Tsukasa Miyajima, Georg 
Hoffstaetter, and Dave Rice. This work is supported 
by the NSF grant PHY-0131508 and NSF/NIH-NIGMS 
award DMR-0225180. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 

m 

[9] 
[10] 

[11] 
[12] 

[13] 
[14] 



J. Billen and L. Young, Los Alamos Laboratory Technical 
Report No. LA-UR-96-1835, 1996 (rev. 2005). [15 
K. Flottmann, Astra - a space charge tracking algorithm, 
available at http://www.desy.de/~mpyflo 
GPT - general particle tracer, Pulser Physics, 
http://www.pulser.nl/gpt [16 
M. Ferrario, J. E. Clendenin, D. T. Palmer, J. B. Rosen- 
zweig, and L. Serafini (2000), SLAC-PUB-8400. 
L. Giannessi and M. Quattromini, Phys. Rev. ST Accel. 
Beams 6, 120101 (2003). [17 
D. P. Rusthoi and K. R. Crandall, Los Alamos Labora- 
tory Technical Report No. LA-UR-97-886 (1997). 
MAFIA, CST GmbH, Darmstadt, Germany, [18 
http : //www. est . de 

D. P. Grote, A. Friedman, and I. Haber, in Proceedings [19 
of the 1996 Computational Accelerator Physics Confer- 
ence (AIP Conference Proceedings 391, Williamsburg, [20 
VA, 1996), p. 51. 

A. Cianchi, D. Alesini, A. Bacci, M. Bellaveglia, R. Boni, 

M. Boscolo, M. Castellano, L. Catani, E. Chiadroni, [21 

S. Cialdi, et al., Phys. Rev. ST Accel. Beams 11, 032801 

(2008). 

V. Miltchev, K. Abrahamyan, G. Asova, J. Bahr, G. Dim- [22 
itrov, H.-J. Grabosch, J. Han, M. Krasilnikov, D. Lipka, 
A. Oppelt, et al., in Proceedings of the 2004 PEL Con- [23 
ference (Trieste, Italy, 2004), p. 399. 
Mini-Workshop On Characterization Of High Bright- 
ness Beams, DESY, Zeuthen, Germany, 2008, [24 



https : //indico .desy . de/ conf erenceDi splay .py?conf Id=806 

B. M. Dunham, C. K. Sinclair, I. V. Bazarov, Y. Li, [25 
X. Liu, and K. W. Smolenski, in Proceedings of the 
2007 Particle Accelerator Conference (Albuquerque, NM, [26 
2007), p. 1224. 

J. Billen and L. Young, Los Alamos Laboratory Technical 
Report No. LA-UR-96-1834 (2000). 

Belomestnykh, V. Shemelin, K. Smolenski, and 
V. Veshcherevich, in Proceedings of the 2007 Particle Ac- 



celerator Conference (Albuquerque, NM, 2007), p. 2331. 
D. G. Ouzounov, I. V. Bazarov, B. M. Dunham, C. K. 
Sinclair, S. Zhou, and F. Wise, in Proceedings of the 
2007 Particle Accelerator Conference (Albuquerque, NM, 
2007), p. 530. 

I. V. Bazarov, D. G. Ouzounov, B. M. Dunham, Y. Li, 
X. Liu, R. E. Meller, J. Sikora, C. K. Sinclair, F. W. 
Wise, and T. Miyajima, Phys. Rev. ST Accel. Beams 

II, 040702 (2008). 

S. G. Anderson, J. B. Rosenzweig, G. P. Le Sage, and 
J. K. Crane, Phys. Rev. ST Accel. Beams 5, 014201 
(2002). 

C. Lejeune and J. Aubert, Adv. Electron. Electron Phys. 
Suppl. A 13 (1980). 

S. Agostinelli et al. (GEANT4), Nucl. Instrum. Meth. 
A506, 250 (2003). 

C. Gulliford, Cornell University, In- 

ternal Report, 2006. Available at 

http: //www. lepp . Cornell . edu/~ib38/reu/06/gullif ord. 
I. V. Bazarov, B. M. Dunham, Y. Li, X. Liu, D. G. 
Ouzounov, C. K. Sinclair, F. Hannon, and T. Miyajima, 
J. Appl. Phys. 103, 054901 (2008). 

L. M. Young and J. Billen (2003), particle Accelerator 
Conference (PAC 03) 12-16 May 2003, Portland, Oregon. 
S. B. van der Geer, O. J. Luiten, M. J. de Loos, 
G. Poplau, and U. van Rienen (2003), DESY-TESLA- 
2003-04. 

M. P. Stockli, R. F. Welton, and R. Keller, Rev. Sci. 
Instrum. 75, 1646 (2004). 

I. V. Bazarov and C. K. Sinclair, Phys. Rev. ST Accel. 
Beams 8, 034202 (2005). 

I. V. Bazarov, B. M. Dunham, F. Hannon, Y. Li, X. Liu, 
T. Miyajima, D. G. Ouzounov, and C. K. Sinclair, in 
Proceedings of the 2007 Particle Accelerator Conference 
(Albuquerque, NM, 2007), p. 1221. 



